Staggered and extreme localization of electron states in fractal space 
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We present exact analytical results revealing the existence of a countable infinity of unusual 
single particle states, which are localized with a multitude of localization lengths in a Vicsek fractal 
network with diamond shaped loops as the 'unit cells'. The family of localized states form clusters 
of increasing size, much in the sense of Aharonov-Bohm cages [J. Vidal et al, Phys. Rev. Lett. 81, 
5888 (1998)], but now without a magnetic field. The length scale at which the localization effect for 
each of these states sets in can be uniquely predicted following a well defined prescription developed 
within the framework of real space renormalization group. The scheme allows an exact evaluation of 
the energy eigenvalue for every such state which is ensured to remain in the spectrum of the system 
even in the thermodynamic limit. In addition, we discuss the existence of a perfectly conducting 
state at the band center of this geometry and the influence of a uniform magnetic field threading 
each elementary plaquette of the lattice on its spectral properties. Of particular interest is the case 
of extreme localization of single particle states when the magnetic flux equals half the fundamental 
flux quantum. 

PACS numbers: 73.20.Fz, 73.22. Dj 



I. INTRODUCTION 

Interplay of lattice topology and quantum interference 
effects is known to give rise to exotic electronic spectrum 
in solid systems that has been studied in details over 
several decades by now. The subject is still being pursued 
with vigor and with an aim to achieve comprehensive 
control over coherent transport in low dimensions. 

In quantum interference and related transport mecha- 
nism a pivotal role is played by the localization of elec- 
tronic eigenstates in presence of disorder. Such localiza- 
tion, known as the Anderson localization [l| upholds a 
central result that, in one dimension with arbitrary dis- 
order, all the single particle states will be exponentially 
localized, and the same was shown to ring true in two 
dimensions as well [|| ■ Since then, extensive research has 
been undertaken to understand the fundamentals of lo- 
calization effects, studies ranging from electronic states 
in random lattice models Q-Hi to the Anderson localiza- 
tion of light I3.fl0j. spin freezing in one dimensional semi- 
conductors , and localization in optical lattices [l2j , to 
name a few. Matter waves can also be localized in deter- 
ministic potentials sharing certain features of random dis- 
order (l3L II4I ] . Recent experiments reveal the Anderson 
localization of non- interacting Bose-Einstein condensates 
in one dimensional matter waveguides, where the random 
potential has been generated by laser speckles [l5| . Simi- 
lar experiments have also been reported to study the An- 
derson localization i n op tical lattices [TH, [TtJ , and in the 
cases of microwaves [18[ and of classical waves in weakly 
disordered one dimensional stack of meta-materials fl9j . 

Variations of the classic Anderson localization are also 
well known by now. Isolated de-localized (extended) sin- 
gle particle states exist, even in a disordered one di- 
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mcnsional chain of atomic potentials, resulting out of 
a kind of spatial correlation (2^, HJ, in 1-d quasiperi- 
odic chains [22l-[2r|. or, in certain kinds of deterministic 
fractal geometries [27|]-|29|]. Crossover from an insulating 
to a metallic spectral behavior in correlated disordered 
two-le gged ladder networks have also been reported re- 
cently^, HI . 

A curious point, apparently gone unnoticed or un- 
appreciated so far is that, while a precise determination 
of the eigenvalues corresponding to the extended single 
particle states is possible in the above cases of correlated 
and deterministic disorder, the task seems to be practi- 
cally impossible when it comes to an exact evaluation of 
eigenvalues of the localized states in a random or even 
a dcterministically disordered system in the thermody- 
namic limit. It should be appreciated that, though a di- 
rect diagonalization of the Hamiltonian for a finite size of 
the system yields eigenvalues of the localized states (for 
a disordered or a deterministically disordered system), 
there is no apriory reason to assume that these eigenval- 
ues remain in the spectrum when the system grows in 
size, and tends to infinity. In fact, for a deterministic 
fractal geometry that offers a singular continuous spec- 
trum, it is almost impossible to hit the exact eigenvalues 
corresponding to the states that will finally be localized 
on an infinite lattice. To the best of our knowledge, this 
issue remains unaddressed so far in the literature. 

Can one really identify the localized states and extract 
the corresponding eigenvalues for a deterministically dis- 
ordered system? In the present communication we ad- 
dress ourselves this question, and take up the task of 
critically exami ning the spectral properties of a Vicsek 
fractal network [32| consisting of diamond shaped loops 
within a tight binding formalism. While looking for the 
localized state eigenvalues and the nature of localization 
are indeed the major factors driving this work, other in- 
terests in such a study are related to the general spectral 
character and magneto-transport in such systems. The 
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motivation behind the latter part of the work may be 
summarised as follows. A diamond- Vicsek network (see 
Fig. [T]) provides an interesting geometry in which the 
'open' character of a typical Vicsek pattern is preserved 
along with the presence of closed loops in shorter scales 
of length. This is in marked contrast to the much stud- 
ied Sierpinski gasket [HI, [H| , which is a closed structure, 
or to the other open tree fractals [35j or even an alterna- 
tive version of the Vicsek fractal without any local closed 
loops (36| . The presence of these loops effectively gener- 
ates a longer ranged interaction between the atomic sites 
occupying the various vertices, and its effect on the elec- 
tron localization or de-localization is worth studying. 

Secondly, linear arrays of diamond networks have al- 
ready drawn considerable attention in recent years in the 
context of charge and spin transport properties (37j-[3^|. 
being shown to behave as a flux controlled n- or p-type 
semiconductor [37|. or as a prospective candidate of an 
elegant spin filter |38l. l39j . The influence of a topological 
variation in the arrangement of such loops on the spec- 
tral properties of the system is thus worth investigating, 
both from the standpoint of fundamental physics, and 
from the perspective of device technology. We choose 
such a deterministic geometry to make an analytical at- 
tack on the system possible. 

We find extremely interesting results in the context of 
localization of electronic states. In the absence of any 
external magnetic field, a countable infinity of localized 
states can be precisely detected with a multitude of local- 
ization lengths. One can work out an exact mathematical 
prescription to specify the length scale at which the onset 
of localization takes place. The localization can in prin- 
ciple, be delayed (staggered) in position space and the 
corresponding energy eigenvalues can be exactly evalu- 
ated following the same prescription based on a real space 
renormalization group (RSRG) method. In addition, it is 
shown that for a given set of parameters, the center of the 
spectrum corresponds to a perfectly extended cigenstate, 
with the parameters of the Hamiltonian exhibiting a fixed 
point behavior. Switching on a magnetic field opens up 
gaps in the spectrum in general, and even leads to an ex- 
treme localization of all the single particle states in the 
sense of formation of the Aharonov-Bohm cages [io[ • 

In what follows we describe the results. In section II, 
the model and the mathematical method of handling the 
problem are presented. Section III and IV include the re- 
sults and their analyses related to the spectral properties 
without and with the magnetic field respectively. The 
two terminal transport study is carried out in section V 
and in section VI we draw our conclusions. 



II. THE SYSTEM AND THE MATHEMATICAL 
FORMULATION 

A. The Hamiltonian 

We refer to Fig. [IJa) which illustrates the second gen- 



eration of a Vicsek geometry with diamond shaped loops. 
As, at one stage we shall be considering the effect of a 
magnetic field on the spectral properties, we show in the 
figure the flux distribution. Each plaquette is threaded 
by a uniform magnetic flux $. The atomic sites are as- 
signed different status depending on their positions and 
neighborhood in the lattice, namely, the sites in the bulk 
are marked as l B\ while 'A' refers to the sites sitting at 
the edges. The magnetic field breaks the time reversal 
symmetry along the edges of every diamond. The Hamil- 
tonian, in the tight binding formalism is written as, 

H = X>|i><ii+x;[ i f > 01 + ^~ i9ii u) m ] a) 

where, q is the on-site potential at the i-th site, and 
has a value 6 a or es depending on whether its an 'edge' 
site or a 'bulk' one. The uniform nearest neighbor hop- 
ping amplitude is tij = t along the edges, and = X 
when i and j refer to the opposite vertices, connected 
by a diagonal. Thus we keep the provision of including 
hopping beyond the nearest neighbors. 6ij is the Pcierls' 
phase [41| given by, 6*y = 27r<E> / 4<i>o f° r hopping along an 
edge. <&o = he/ e is the fundamental flux quantum. From 
symmetry considerations, Oij = when the electron hops 
across a diagonal, that is, when tij = X. 

B. The RSRG scheme 

An elegant way of handling such self-similar systems 
is to use the real space renormalization group (RSRG) 
method [42j where one can decimate out a subset of 
atomic sites from the original lattice to get a scaled ver- 
sion of it (Fig.^b)). This is easily done by writing down 
in details the set of difference equations, 

(E-a) ^ = 4>j ( 2 ) 

3 

where, ipi denotes the amplitude of the wave function at 
the i-th site, and 0^ is the Pcierls' phase associated with 
the hopping matrix element connecting the i-th and the 
j-th sites. We begin with nearest neighbor hopping only 
(that is, we set A = at the beginning). However, such a 
decimation automatically generates the second neighbor 
hopping across the diagonals of an inflated diamond as 
shown by the dotted line in Fig. [IJb). The range of inter- 
actions of course does not increase beyond this on further 
renormalization. It is therefore advisable to retain A in 
the Hamiltonian from the very beginning. One can eas- 
ily compare the results obtained by switching A on or off. 
The recursion relations of the on-site potentials and the 
hopping integrals are provided below. 

('a = e A + [ptf + p*h + a\i] 

('b = e B + 2 [ptf + p*t b + aXi] 

t'f = /3Ai; t' b = 0*X i; A' = 7 Ai (3) 
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FIG. 1: (Color online) (a) Schematic view of the second generation of an infinite diamond- Vicsek network with each diamond 
plaquette threaded by a uniform magnetic flux The edge sites are named l A' (blue circles) and the bulk sites are marked as 
l B' (orange circles), and the arrow indicates the direction of the forward hopping. We exclude any second neighbor hopping at 
the beginning, (b) Renormalized version of (a) with the dotted lines indicating the diagonal hopping which is generated due 
to renormalization. 



where, a = [(E-e B )\i]/&, P = [(E-e B )tf+X 2 t h }/&, 
(3* = \{E - l B )t b + \ 2 i f }/&, 7 = A 1 A 2 /^ 3 ; 6 = {E - 

with, e B = e B + w*tj + wt b - if = utj + vt b ; i b = 
vtf + ut b ; X2 = X + wtf + w*t b . 

Here, u = \{E - e s )Ai]/&; v = AA1/6; w = \{E - 
i B )t f + At 6 ]/&; w* = [(E - e B )h + At/]/&; 6 = (E - 

with, e B = e B + ptf +p*t b ; Ai = A +p*tf + pt b 
andp=[(E-e A )t b + Xt f ]/^;p* = [(E -e A )t f + Xt b }/^; 
^^(E-e A ) 2 -X 2 . 

In the above expression, tf = t b = te l6 , where 8 = 
27r<I>/4$o is the constant Pcierls' phase. The above recur- 
sion relations are then used to obtain information about 
the local density of states (LDOS) at specific sites of the 
system, and the character of the single particle states, as 
discussed below. 



III. SPECTRAL PROPERTIES WITH ZERO 
MAGNETIC FIELD 



A. Local density of states in zero magnetic field 
and with A = 



Using the standard decimation procedure [42(, the 
LDOS at the edge (A) and the bulk (B) sites can eas- 
ily be obtained through the local Green's functions. For 
simplicity we present in Fig. [UJa) the LDOS at a £>-site 



only, given by, 



P (B) {E) = lira 



7)-S-0 



--Im{G (s) (£; + ^)} 



(4) 



where, G^ B \E+ivi) = (E+irj — eg)^ 1 , e* B being the fixed 
point value of the relevant on-site potential at the i?-sitc, 
obtained by iterating Eq. ([3]). We have set e^i = £5 = 0, 
t = 1 and A = 0, and there is no magnetic field (i.e., 
$ = 0). The LDOS shows a dense packing of eigenstates 
over a finite range of energy centered at E = 0. We have 
minutely examined this continuum by fine scanning an 
energy interval around E = 0, and show it in Fig. [UJb). 
The continuum seems to persist. In the neighborhood 
of the band-center and within the apparent 'continuum', 
the hopping integral remains non-zero over a substantial 
number of RSRG iteration steps. The number of such 
steps n depends on the chosen energy and indicates that 
the corresponding eigenfunction is either of an extended 
character, or, at least, has very large localization length. 
This aspect will be further discussed in the following sub- 
section in relation to the so called staggered localization 
effect. 

A particularly interesting state is the band cen- 
ter, viz., E = 0, where the entire parameter space 
{sa, e_B, A} exhibits a one cycle fixed point behav- 
ior, viz., {e A (n + l),e B (n + l),i(n + l),A(n + 1)} = 
{eA(n), e B (n), t(n), A(n)} for n > 1. n stands for the 
RSRG iteration step. A at this special energy remains 
zero throughout the iteration. We conclude that the 
eigenstate at the band-center is definitely extended, but 
is of a non-Bloch character. The general behavior of the 
hopping integrals under successive RSRG iterations is 
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FIG. 2: (Color online) (a) LDOS-S plot at the bulk sites 
(B-type) of an infinite diamond- Vicsek network in absence 
of the magnetic field ($ = 0). The other parameters are 
€a — £b — 0, t = 1 and A = 0. (b) A highlighted version of 
(a) around the center E = 0. 



suggestive of the fact that this central extended eigen- 
state might be flanked on either side by a countable in- 
finity of eigenstates which belong either to the extended 
category, or have very large localization lengths. 



B. Exact construction of eigenstates 

The inherent self-similarity of the deterministic fractals 
allows for the construction of exact distribution of ampli- 
tudes of the eigenstates, by suitably exploiting Eq. 
Previous attempts in this regard have unfolded extended 
non-Bloch states ( atypically extended states) in the cases 
of an open loop Vicsek fractal 1431) or a closed loop di- 
amond hierarchical geometry [44|. The present lattice 
offers a richer spectrum, allowing one to explicitly con- 
struct localized states extending over clusters of lattice 
points of various sizes on the parent lattice. The planar 
extent of such clusters depends on the eigenvalue cor- 
responding to the localized state, and can be small or 
enormous. 

To elaborate, let us consider the solutions of the equa- 
tion, 

E = e B (n) - 2A(n) (5) 
where, n refers to the stage of rcnormalization. This is 



in general, a polynomial equation in E. The zero's of the 
polynomial will be eigenvalues of the infinite system if, 
and only if, with them one can satisfy Eq. ([2]) locally at 
every vertex of the lattice, even when the lattice grows 
infinitely large. This task can be accomplished by trying 
to draw a non-trivial distribution of amplitudes for an en- 
ergy obtained from Eq. ([5]) on the undecimated vertices 
of an n-step renormalizcd lattice, and then trying to fig- 
ure out the amplitude distribution on the original lattice 
at the bare length scale. Let us discuss two specific cases 
at first. 

Case I: We begin with the un-renormalized lattice. 
Now n = 0, and with €a = (b = and A = the solution 
of Eq. ([5]) is E = 0. One can construct an eigenfunction 
for E = with amplitudes equal to ±1 distributed alter- 
nately at the '-B' sites along the major X- and Y-axes. 
The difference equation, viz., Eq. ([2]) can then easily be 
satisfied for all other intermediate vertices using the val- 
ues or ±1. 

Case II: The above idea can indeed be extended to 
higher values of n, as we demonstrate in Fig. [3Ja) for 
n = 1, and discuss below. Let us extract the roots of 
the Eq. ([5]) for n = 1. The roots are E = and ±\/6 
for €a = £_b = 0, t = 1 and A = initially. We explain 
the construction of amplitudes for E = \/6. The trick in 
this case is to place the values ±1 periodically along the 
major A-and Y-axes, and to assign an amplitude equal 
to zero at every 'edge' (A-type site) on a one step renor- 
malizcd lattice. The amplitudes at the intermediate sites 
of the original lattice are then systematically evaluated 
using Eq. We show it in Fig. [3Ja) on a small por- 
tion of an infinite lattice. The bigger square boxes with 
thick, highlighted edges represent the one step renormal- 
ized lattice. The construction depicted on a smaller scale 
can be extended to see that the distribution holds on a 
lattice of larger spatial extent (Fig.[3Jb)). In fact it holds 
even on a lattice of an arbitrarily large size, where the 
end sites are not actually visible. As, according to our 
earlier argument, we are able to satisfy Eq. ^ locally at 
every vertex while drawing this distribution, E = y/6 is 
definitely an eigenvalue of the infinite system, a fact that 
has been cross-checked by evaluation the LDOS at the A- 
and the -B-sites at this special energy. We get a stable, 
finite value of the LDOS which supports our argument 
above. 

Looking back at Fig. [31 the solid black dots represent 
non-zero values of the amplitude while an open circle 
represents an amplitude equal to zero. Non-zero ampli- 
tudes, represented by solid circles, have- the values equal 
to ±1, ±V6/4 (depicted by the letter a) and ±1/2 (the 
letter b), distributed suitably so as satisfy Eq. © consis- 
tently at every intermediate vertex on the original lattice. 
The grey shaded clusters in Fig. |3ja) embrace the non- 
zero amplitudes only, while every yellow shaded zone is 
surrounded by vertices where amplitudes are zero. The 
significant observation is that, clusters of non- zero am- 
plitude span over a finite distance, but ultimately get 
decoupled from each other on a larger scale of length. 
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FIG. 3: (Color online) (a) Distribution of amplitudes of the wave functions at E = \/6 (obtained by solving Eq. ([EJ for n — 1) 
on a 2nd generation lattice. The dark shaded plaquettes (grey ones) embrace atomic sites with non-zero amplitudes (black 
solid circles) and the light shaded plaquettes (yellow ones) are surrounded by atomic sites with zero amplitudes (white empty 
circles). The initial parameters are chosen as $ = 0, £a = es = 0, t = 1 and A = 0. In the figure, a = -v/6/4 and b = 1/2. 
(b) Distribution of amplitudes of the wave functions at E — v6 on a 3rd generation lattice, the other parameters and symbols 
are same as (a). In the figures (a) and (b), the highlighted thick (blue) lines represent one-step and two-step renormalized 
lattice respectively. 



This can be appreciated if we look at Fig. [3^b) which is 
a larger version of the previous figure. The grey shaded 
clusters are distributed along the principal X- and Y- 
axes, but are decoupled from each other beyond a certain 
extent by the unfilled white boxes. The yellow clusters 
representing amplitude- voids are now seen to span larger 
spatial distances. A similar construction is possible for 
E = — \/6 which is another solution of Eq. ([5]) for n = 1. 



C. Staggered localization 

It is appearent from the above discussion that the 
cigenfunction corresponding to E = ±^/6 will be lo- 
calized in the fractal space, as the spanning clusters of 
different sizes and embracing non-zero amplitudes ulti- 
mately get decoupled from one another. This is easily 
re-confirmed by studying the evolution of the hopping 
integrals under successive RSRG steps. The hopping in- 
tegrals t and A (zero initially, but grows later) remain 
non-zero at the first stage of RSRG (that is, n = 1), in- 
dicating that the nearest neighboring sites on a one step 
renormalized lattice will have a non-zero overlap of the 
eigenfunctions. They start decaying for n > 2 with the 
decay in A(n) taking place at a much slower rate com- 
pared to t(n). This indicates that over larger scale of 
length the corresponding states are ultimately localized, 
but the effect is a weak one. 



This observation immediately leads to an innovative 
way of exactly determining a localized cigenstate on such 
a deterministic geometry. It should be appreciated that 
though it is not unnatural that most of the single par- 
ticle states will be Anderson-localized in the absence of 
any translational order, nevertheless an exact prescrip- 
tion of the determination of any localized eigenvalue is 
not easy to obtain, and has not been reported so far in 
the literature, to the best of our knowledge. We do it 
using the following method. 

We can solve Eq. ([5J in principle, for any n. For ex- 
ample, we have done it explicitly for n = 1, n = 2, n = 3 
and n = 4. With the same set of parameters as dis- 
cussed above, the roots of Eq. ([5]) for n = 2 are, E = 0, 
±VE, ±2.11619, ±0.77508, ±2.98681. As we observe, 
the roots for the n = 1 stage, viz., E = and ±v / 6 
are included in this set for n = 2. E = corresponds 
to the extended state and E = ±v / 6 provide two local- 
ized states we already know. For each of the additional 
roots, viz., E = ±2.11619, ±0.77508, ±2.98681, the hop- 
ping integrals t and A remain non-zero (with considerable 
magnitude) up to the second stage of iteration (n = 2), 
and starts to lose their 'strengths' as the renormalization 
progresses. Finally, for large 'n' the hopping integrals 
become zero. 

The above observation implies that, using a subset 
of energy values extracted at the stage n = 2 (E = 
±2.11619,±0.77508 and ±2.98681), we can work out 
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cigcnfunctions which will span bigger clusters of lattice 
points on the original lattice compared to those obtained 
from n = 1. The states will appear to be 'extended' when 
viewed on a finite diamond- Vicsek fractal at the second 
generation, but will eventually be localized on a lattice 
in the thermodynamic limit. In Fig. 0] we show the dis- 
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cigenstate determined so far resides. The clustering of 
the eigenvalues is shown in Fig. [5] This dense filling of the 




FIG. 5: (Color online) Distribution of energy eigenvalues E 
obtained from the Eq. (0 for different RSRG stage n. The 
central dot (red one) at E = represents the eigenvalue for 
the extended eigenstate. 

eigenvalue spectrum around the center is also reflected in 
the apparent continuum observed in the density of states 

(Fig- EDO). 



FIG. 4: (Color online) Distribution of amplitudes of the 

wave functions at E = 0.77508 (obtained by solving Eq. © IV. SPECTRAL PROPERTIES WITH 

for n = 2) on a 3rd generation lattice. The other parameters NON-ZERO MAGNETIC FIELD 

and symbols are same as FIG. [3] The thick, highlighting blue 

lines represent a two-step renormalized lattice. A. The energy eigenvalue spectrum 



tribution of amplitudes for E = 0.77508, a value that is 
obtained from Eq. ([S]) for n — 2. The enlargement in 
the cluster-size having non-zero values of the amplitude 
in comparison to the n = 1 case (Fig. E^b)) is obvious. 
The spanning clusters finally get decoupled from each 
other, just as it was for the n = 1 case. But, now this 
decoupling occurs at a larger length scale. 

It is now easy to foresee what is going to happen for 
n = 3, 4 and beyond. For any n = £ we will be getting 
roots of Eq. ([5]), subsets of which are solutions of Eq. §5§ 
for n = 1, 2, . . . , £ — 1. For these subsets, the decay 
in the hopping integrals will begin at n > 1, n>2, 
n > £ — 1. For the roots in addition to these, the hop- 
ping integrals lose their strengths and finally decay, from 
n > £. Thus, the latter eigenvalues will correspond to lo- 
calized eigenstates, the localization being delayed (stag- 
gered) in space with localization lengths much larger than 
the previous ones. When mapped back on to the original 
lattice, the amplitudes for these additional roots will be 
found to span clusters of increasing size. The exact size 
of the spanning clusters will be determined by the value 
of n. 

The roots of Eq. §5§ are found to cluster around the 
value E = symmetrically, and tend to densely fill the 
neighborhood of E = 0, at which the single extended 



We have obtained the energy eigenvalue distribution 
(Fig. D as a function of the magnetic flux $ enclosed in 
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FIG. 6: (Color online) Energy eigenvalue spectrum of an 
infinite diamond- Vicsek fractal as a function the magnetic 
flux $. We have chosen 6a — es = 0, t = 1 and A = 0. 

each basic plaquette for an infinite size diamond- Vicsek 
fractal. We have examined the formation of the bands 
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and the gaps with the variation of magnetic flux i>. To 
obtain the energy spectrum, we have calculated the local 
density of states (LDOS) at both 'A' and l B' sites by 
fixing the value of the energy E and varying the mag- 
netic flux $ from —1 to 1, repeated the above process for 
different values of energy E, and picked up those values 
of energy E and magnetic flux <£> for which we get a non- 
zero LDOS cither at an A-site or a £>-site. Thus Fig. [5] 
is representative of an infinite lattice. 

In Fig. [IJJ we can clearly observe the formation of mul- 
tiple bands and gaps, and how a variation of magnetic 
flux $ leads to band overlapping. The band crossing is 
maximum at the center (around $ = 0) and the density 
of allowed energy eigenvalues is large is this area. As we 
shift form $ = on either side, there is thinning of al- 
lowed energy eigenvalues. And finally at $ = <£>o/2, only 
four energy eigenvalues are allowed indicating an extreme 
localization of the electronic states which is discussed in 
details in the next subsection. Fig. [6] corroborates this 
last observation, though due to limit of resolution of the 
diagram, the four eigenvalues are not clearly seen there. 
This is resolved in the next diagram (Fig. [3b)). 



B. Extreme localization of the electronic states 

In absence of magnetic field, there was clearly a non- 
zero value of LDOS at the center of the spectrum (around 
E = 0) (Fig. Ufa)). As soon as the magnetic field is 
switched on, a wide gap opens up in the LDOS spectrum 
around E = (Fig. [TJa)). The gap becomes wider as we 
increase the value of magnetic flux $, finally leading to 
an extreme localization of electronic states (Fig. 0(b)) at 
the half flux quantum (i.e., <f> = $o/2). The four separate 
lines in Fig. (at E = ±y/2 and at E = ±2) corre- 
spond to four highly degenerate localized states pinned 
at the 'A'- and 'B'-sites respectively. This observation 
is in accordance with Vidal et al.'s result [4(|, where the 
Aharonov-Bohm caging of the localized orbitals under 
the action of an external magnetic field was discussed. 

The origin of the four localized state eigenvalues above 
can be easily explained if one appreciates that at $ = 
$o/2, the effective coupling between the sites at the ver- 
tices of an elementary rhombus, viz., t e ^ = [2t 2 /(E — 
ca)} cos(7t<I>/<I>o), becomes equal to zero. In that case one 
is left with only two types of atomic sites, decoupled from 
each other, and having effective on-site potential energies 
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and t = 1, the localized states are obtained by 

eff 
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setting E = e e £ and E = e'jf, which yield the values 
E 



±y/2 and ±2 respectively. These are the energy 
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FIG. 7: (Color online) (a) LDOS-B plot at the bulk sites (B- 
type) of an infinite diamond- Vicsek network with $ = <3>o/4, 
and (b) is for <E> = $o/2. The red lines correspond to LDOS 
at B-sites and the green lines correspond to LDOS at ^4-sites. 
We have chosen ca = es = 0, t = 1 and A = 0. 



eigenvalues at which extreme localization is observed, as 
shown in Fig. 0L). 



V. TWO TERMINAL CONDUCTANCE FOR A 
FINITE LATTICE 

To get the two terminal conductance for a finite size 
diamond- Vicsek fractal, we attach the system between 
two semi-infinite one-dimensional ordered metallic leads, 
namely, the source and the drain. The leads, in the tight 
binding model, are described by a constant on-site po- 
tential e/ and a nearest neighbor hopping integral t\. We 
then successively rcnormalize the system to reduce it into 
an effective di-atomic system [43^ , consisting of two renor- 
malized atoms, each having an effective on-site potential 
equal to U and with an effective hopping integral V be- 
tween them. The transmission coefficient across the ef- 
fective dimer is given by the well known formula [45[ , 
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T = 



4 sin ka 

[(Mu - M 2 i) + (Mn - M 22 )cosfca] 2 + [(M n + M 22 )sinfca] 2 



(7) 



r 



where, Mn 



-Mi 2, M 2 



Vtj 

_U 

V 



V 

-, M 12 



, M 21 



are the matrix elements of the 



transfer matrix for the effective di-atomic system, and 
cosfca = (E — ei)/2ti, 'a' being the lattice constant and 
is taken to be equal to unity throughout the calculation. 

In Fig. [HI we have shown the two-terminal transmission 
characteristics for a 3rd generation system for different 
values of magnetic flux For $ = 0, the system exhibits 
a continuous high transmission window over a region at 
the center (Fig. Efa)). This is due to the fact that on ei- 



localization lengths extending much beyond the 3rd gen- 
eration fractal. 

As we tune the magnetic flux to a non-zero value, e.g., 
$ = $o/4, the transmittivity of the system drastically 
decreases (Fig. [5th)) and with increase in the value of 
the transmission coefficient T decreases more and more 
and finally at $ = $o/2, the system becomes completely 
opaque to an incoming electron. So by fixing the Fermi 
level of the electron to a particular energy, say at E = 0, 
one can easily transform the system from a conducting 
one to an insulating one by tuning the external magnetic 
flux $ suitably. 
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VI. CONCLUDING REMARKS 

In conclusion, we have examined the energy spectrum 
of a Vicsck geometry consisting of diamond shaped loops. 
The major result is that we have been able to identify 
a countable infinity of localized eigenstates displaying a 
multitude of localization lengths. A prescription is given 
for an exact determination of the eigenvalues correspond- 
ing to all such states, a problem that is far from trivial 
in the case of a deterministically disordered system. The 
localized states span across the fractal space in clusters of 
increasing sizes, the size being precisely controlled by the 
length scale at which the energy eigenvalue is extracted. 
The onset of localization can be exactly predicted from 
the stage of RSRG. In addition, the application of a uni- 
form magnetic field perpendicular to the plane of the 
fractal is found to produce gaps in the energy spectrum. 
A special value of the magnetic flux viz., <f> = $o/2 is 
shown to lead to an extreme localization of the electron 
states as well. The results are corroborated by the den- 
sity of states calculations, and the valuation two terminal 
magneto-transport . 



FIG. 8: (Color online) Transmission characteristics for a 3rd 
generation system, (a) is for $ = and (b) is for $ = $o/4. 
The other parameters are €a ~ ~ 0, t — 1 and A = 0. The 
lead parameters are ei — and ti = 1. 

ther side of E = 0, the energy eigenvalues become quite 
densely packed. The corresponding cigenfunctions have 
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